home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 001 / pn.arc / PNT1.DAT < prev    next >
Encoding:
Text File  |  1985-08-12  |  1.8 KB  |  59 lines
  1. "            x                  PN(x)           abs(error)
  2. "
  3.        -12.73      .0
  4.        -12.72      .0
  5. "
  6.     -9.01327  1E-19
  7.     -7.34880  1E-13
  8. "
  9.         -6.70602  1E-11
  10.         -5.99781  1E-09
  11.         -5.19934  1E-07
  12.     -4.75342  1E-06
  13.         -4.26489  1E-05
  14.         -3.71902  1E-04
  15. "
  16. "       The values computed for PN(x),  with x very small,  are
  17. "       compared with tabulated values of PN(-x) = QN(x) in The
  18. "       Handbook of  Mathematical Functions  [HMF: Table 26.6].
  19. "       Errors within 1E-5 are usually quite acceptable.   Note
  20. "       the very high relative errors out beyond -6.0, however.
  21. /            z                  PN(x)           abs(error)
  22. "
  23.     -0.14         .4443299952
  24.     -0.02         .4920216863
  25.     -0.00251     .499
  26. "
  27.     -3.647656264E-11 .499999999999
  28.     -3.647514842E-11 .5    {-2.579182479E-11 * sqrt2}
  29.      0.0             .5
  30.      3.647514842E-11 .5
  31.      3.647656264E-11 .500000000001
  32. "
  33.      0.00251     .501
  34.      0.02         .5079783137169
  35.      0.14         .5556700048059
  36. "
  37. "       Midrange values, near .5, check on the smoothness of
  38. "       the crossover from negative to positive values of x.
  39. "       Test values are from [HMF:Table 26.1] and ERFT1.PAS.
  40. /            x                  PN(x)           abs(error)
  41. "
  42.     2.32635     .99
  43.     3.09023     .999       
  44.         3.71902     .9999
  45.     4.26489        .99999
  46.     4.75342     .999999
  47.     5.19934     .9999999
  48.     5.61200     .99999999
  49.         5.99781     .999999999
  50.         6.36134     .9999999999
  51.         6.70602     .99999999999
  52.         7.03448     .999999999999
  53. "
  54. "       Near the tail, as PN(x) approaches 1.0,  values of 1-QN(x)
  55. "       from [HMF: Table 26.6] are used to see how well PN(x) does
  56. "       despite adjustments by 1.0 within the PN procedure.
  57. "
  58. "PNT1> End of test.
  59.